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Abstract 

The experiment aiming at the simultaneous determination of the two trans- 
versal polarisation components of electrons emitted in the decay of free, po- 
larised neutrons is in progress at the Paul Scherrer Institute (Villigen, Switzer- 
land). The non-zero value of R coefficient, proportional to the polarisation 
component, which is perpendicular to the plane spanned by the spin of the 
decaying neutron and the electron momentum, would prove a violation of time 
reversal symmetry and thus physics beyond the Standard Model. The planned 
accuracy of the measurement is of order 0.005. To reach this value, the system- 
atic effects in the experiment have to be controlled on a similar level of accuracy. 
The emphasis of this master's thesis is put on the search of systematic effects 
by the means of dedicated Monte Carlo simulation, based on extended Geant4 
package. Implementation details are discussed and the new added features are 
tested. Finally, the (3 decay asymmetry induced systematic effect, resulting in 
false contribution to .R-coefficient is recognised and investigated. 



Streszczenie 

W Instytucie Paula Scherrera (Villigen, Szwajcaria) prowadzony jest ekspery- 
ment majacy na celu jednoczesny pomiar obu poprzecznych skladowych po- 
laryzacji elektronow emitowanych w rozpadzie swobodnych, spolaryzowanych 
neutronow. Niezerowa wartosc wspolczynnika R, proporcjonalnego do tej skla- 
dowej polaryzacji, ktora jest prostopadla do plaszczyzny tworzonej przez spin 
neutronu i ped elektronu, swiadczylaby o istnieniu procesu lamiacego symetrie 
wzgledem odwrocenia czasu, wykraczajqcego poza obecnie uznawan^ strukture 
Modelu Standardowego. Aby moc osiagnac planowana. dokladnosc pomiaru 
(0.005), konieczne jest oszacowanie mozliwych efektow systematycznych. W 
ponizszej pracy magisterskiej nacisk polozono na poszukiwanie efektow sys- 
tematycznych przy pomocy specjalnie do tego celu stworzonej symulacji Monte 
Carlo, opartej na rozszerzonym pakiecie Geant4. Omowione sa szczegoly im- 
plementacji oraz wyniki testow nowych, dodanych do pakietu funkcji. Opisano 
tez i przeanalizowano odkryty efekt systematyczny, spowodowany asymetriq 
rozpadu (3 i wnoszacy falszywy wklad do mierzonego wspolczynnika R. 
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Chapter 1 



Introduction 



From the nuclear physicist's point of view next few years are going to be really 
interesting. On the one hand, the Standard Model (SM) has still an impressive 
predictive power and so far in particle physics there has been a total agreement with 
experimental evidence. On the other hand, SM is not believed to be a complete 
theory. Nowadays there are plenty of much more elegant new theories or at least 
extensions of the SM. Moreover, there is one huge discrepancy between its predictions 
and astronomical observations, namely the baryon asymmetry of the Universe which 
is by a factor 10 8 larger than expected. In other words, according to the SM, large 
concentrations of baryon matter, such as galaxies, could not have even existed and 
the Universe should have been more like a uniform "sea" of radiation. 

One of very few reasonable explanations of this paradox requires the breaking 
of combined charge conjugation and parity symmetry (CV). Assuming that CVT 
symmetry (combined time-reversal symmetry T and CV) is conserved, what is re- 
quired by all renormalisable quantum field theories, CV violation implies T violation. 
Both CV or T-violating processes have been so far observed only in the neutral K 
and B meson systems and their mechanism has been already included in the SM by 
introducing a quark mixing mechanism. However, to explain the problem of baryon 
asymmetry additional sources of CV or T breaking have to be found. 

Measurements of vector correlations in particle decays and searching for electric 
dipole moments of particles belong to the most promising ways to discover new T- 
violating processes not predicted by the SM. For both types of experiments, physics 
of cold neutrons is especially interesting, due to the availability of high intensity 
polarised beams. The experiment aimed at the determination of the i?-correlation 
parameter (mixed product of neutron spin, electron momentum and electron spin) 
in the /3~ decay of free, polarised cold neutrons will start data taking this summer 
at the Paul Scherrer Institute (Villigen, Switzerland). It is going to be the first such 
a measurement for the decay of free neutrons, bearing a potential to detect either 
a non-standard value, inconsistent with zero or to provide important constraints for 
the T-violating scalar and tensor couplings in the semileptonic weak interactions. 

The subject of this thesis is the analysis of experimental uncertainties in the R- 
correlation experiment. The measurement is planned to achieve accuracy of 0.005, 
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with the main contribution to the error coming from the counting statistics. In order 
to perform detailed study of possible systematic effects, a dedicated Monte Carlo 
simulation (based on Geant4 package) has been created and exploited. 

Before the main goal could have been completed, the following intermediate steps 
had to be fulfilled: 

• Unification of all existing parts of the source code in one general simulation of 
the whole experiment 

• Implementation of new geometry of the experimental setup 

• Modification of Geant4 libraries to include electron polarisation transport and 
the spin dependent neutron f3~ decay 

• Further software development and optimisation 

• Actual simulations and data analysis 

The next chapter briefly describes the theoretical aspects of the neutron f3~ decay, 
which are essential for understanding the principle of the experiment. In addition, 
interactions of electrons with matter and foundations of the polarisation theory are 
sketched, especially the formulae used in the simulation. The third chapter contains 
detailed information about the measurement itself and the description of the experi- 
mental apparatus. In further chapters, the Monte Carlo simulation and its results are 
presented and discussed. Appendices contain practical information on the software 
usage. 



Chapter 2 
Theory 



2.1 The neutron j3 decay and angular correlations 

The neutron is about 0.2% more massive than a proton, which translates to an energy 
difference of 1.29 MeV. Therefore, from the energy conservation law, it is possible for 
a free neutron to decay with the emission of an electron and an electron anti-neutrino 



n 



p + e~ + u e Q = 782.2 ±0.1 keV. 



Because it is a three-particle decay, produced electrons have continuous distribution 
of momentum and energy (Fig. I2.1JI given by 



Wi(p e ) oc F(Z,p e )p e (p max -p e ) , 

2-jia 1 
W 2 (E e ) oc — ^p e E e {E mc 

P l-e~~ 



E P 



(2.1) 
(2.2) 



where F(Z,p e ) is the Fermi function, /3 = ^-, E n 
and a is the fine structure constant. 
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Figure 2.1: The observed momentum and energy distributions for the (3 electron 
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The decay of neutrons involves the weak interaction (see Fig. I2.2|) . which, accord- 
ing to the theoretical description embedded in the SM (Feynman and Gell-Mann, 
1958), has the strict vector-axial form (V — A). To find T- violating effects in f3 




e 




W 

Neutron decay 
via weak 
interaction 



Figure 2.2: The neutron decay as the transformation of one of the neutron's down 
quarks into an up quark. 

decay, one has to measure at least three vector or axial-vector quantities, namely 
momenta and spins of the different particles involved in the decay (e, u e , p, n); the 
quantities, which are all reversed under the time reversal operation. Rotational in- 
variance requires the observables to be scalars or pseudo-scalars, therefore the lowest 
order T- violating combination of spins and momenta appears in the form of the mixed 
triple product. From the experimentally accessible quantities, four such products can 
be formed (see the Tab. I2.1J1 . 



Correlation 


Broken symmetry 


Definition 
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r,v 


J- (p e x a) 


D 
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J' (Pe X Pu e ) 


V 


T,V 


J- (P x a) 


L 


T 


P- (p e x a) 



Table 2.1: T- violating triple products. J is the neutron spin, p e and a are the 
momentum and spin of the electron, respectively, P denotes the proton momentum 
and p Ue stands for the momentum of the neutrino. 

From now on let's concentrate on the planned experiment, in which the momentum 
and spin of electrons from decaying oriented neutrons will be measured. The decay 
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probability distribution suitable for this case has been derived by Jackson [S| 

,(J)- (p e x a) 
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(2.3) 



where 

b 

A = 
G = 

Q = 

N = 
A 

Using A 



, A 2 +A 

' 1+3A 2 



is the Fierz interference coefficient, 
is the decay asymmetry parameter, 



-1, 

o 1+A 

Z 3+A 2 ' 
Xrn e 



E e i 



denotes the coupling constants ratio Ca/Cv- 
-1.267 (see [20]), one retrieves values: 

A = -0.1162, Q = -0.1160. 



The coefficient b vanishes if there are no scalar and tensor couplings, therefore the 
zero value has been used. The formula 12.31 is the consequence of the interaction 
Hamiltonian density for weak currents, as given by Yang and Lee (Zj. As can be seen 
on the Fig. 12 .3^ R and iV parameters are proportional to orthogonal components of 
the transversal electron polarisation. However, only the i?-correlation reveals direct 




Figure 2.3: Directions of spins of the neutron and the electron from the /3~ decay. In 
the .R-correlation experiment the transversal component of the electron polarisation 
will be measured 

sensitivity to the existence of exotic T-violating scalar and tensor couplings in the 
semileptonic weak interactions. 



2.2 Polarisation theory 

The following brief review of the polarisation theory formalism is based on the Refs. 
[H], [IE] and [Hil- The Stokes parameters are defined below for the general case and for 
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the specific case of electrons. When used as a four- vector, the Stokes parameters allow 
the ordinary polarisation-sensitive cross sections to be written in matrix form, which 
is a very convenient representation for the description of polarisation phenomena. 

2.2.1 Stokes parameters 

In quantum mechanics the wave function describing a pure state of polarisation can 
be expanded in a complete set of orthonormal eigenfunctions. For particles of spin ~ 
this expansion contains only two terms, 

Ip = 01*01 + «202- (2.4) 

The wave functions describing pure states may be chosen in the form 
Thus, the general wave function describing the beam is given by 

*=(£)■ (2 - 6) 

In the specific case of the electron, the explicit form of this function can be found by 
solving the Dirac equation. As a result, one can achieve two independent solutions 
(spin "+" and spin "— ") for the electron and two for the positron. Choosing the 
solutions for the electron, one immediately obtains wave functions (|2.5jl . 

In the rest frame of the electron the four-component solutions of the Dirac equa- 
tion are reduced to two-component spinors, leading to the following expressions for 
the expectation value of the unit matrix and the expectation values of the Pauli spin 
operators: 

I = ty | I | ip) = ( a \ a* 2 ) I Q 1 I f ^ ) = a x a\ - a 2 a* 2 , 
Pi = (ip \ a z \ ip) = ( a{ a 2 ) ( o -l ) ( ^ ) = aia * ~ a2a 2 ' 
P2 = (0 I cr x I 0) = ( a* a* 2 ) ( , n ) ( * I = aia* 2 + Oa<A , 



1 J \ a 2 

Ps = (0 I °y I 0) = ( Oj a *2 ) f j ) ( 02 ) = '^ aia * 2 ~ a2a *" > ' 

This is the set of four so-called "Stokes parameters" which represent observables and 
completely characterise the electron in a pure polarisation state. As it can be seen 
from the definition, the physical interpretation of the parameter I is quite obvious, 
it is simply the total beam intensity. Moreover, if the electron is transformed to its 
rest frame, P = (Pi,P 2 ,P3) can be considered as the spin direction. Note, that the 
beam in a pure state is completely polarised, therefore \P\ = 1. 
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2.2.2 Mixed states 

A partially polarised beam of electrons cannot be represented by a single wave func- 
tion, but by an incoherent "ensemble" of pure states, each characterised by its own 
wave function. In order to describe such a case, the density matrix formalism is used 
(for detailed introduction see |Hj). The density matrix p is a 2 x 2 hermitian matrix 
with positive or zero eigenvalues and trace equal 1. For the special case of a totally 
polarised beam, corresponding to the function (J2.6J) . the density matrix is 

' * 0*02 a* 2 a2 

which can be always brought into the simple form 

1 
= I 

by a unitary transformation. Although, in the most general case, p can be always 
diagonalised, for mixed states both diagonal elements remain nonzero. The resulting 
matrix can be considered as the incoherent superposition of the unpolarised and 
totally polarised beam 

where P (0 < P = p a — pb < 1) is called the degree of polarisation. 

The density matrix can be also expressed using the Stokes parameters 

p= l -{I + p.a) 1 

it is worth mentioning that now, for mixed states, the polarisation vector norm 
\P\= P < 1. 

2.2.3 Matrix representation 

The Stokes parameters are very often written in the form of a four- vector: 



/ I \ 
Pi 
Pi 



I 
P 



The interpretation and some typical examples are given in the table 12.21 and below: 

° represents an unpolarised beam, 
W 



ID 



(^n/\ritLtn a. mnjuni 



Stokes parameter 


Interpretation 


I 
Pi 
P2 
Pz 


Intensity 
|1J Longitudinal polarisation (spin in z direction) 
+i/ / *~ 1 Transversal polarisation (spin in x direction) 
_1 < > +1 Transversal polarisation (spin in y direction) 



Table 2.2: Interpretation of components of the Stokes vector. 



( l \ 




( 1 \ 


±1 




or 




±1 


V ) 




V / 



describe transversal polarisation in z and x directions. 



Because the Stokes parameters are dependent on the polarisation basis one has 
chosen, there exists a transformation matrix M, which relates a Stokes vector in one 
reference frame to the same one in another coordinate system: 



V 

F" 



M 



I 
P 



M 



(2.7) 



For instance, if the second coordinate system is rotated about the direction of propa- 
gation at an angle 9 to the right, then the matrix M is a simple rotation matrix given 
by 

/ 1 \ 

10 

cos# sin# 
\ - sin 9 cos 9 J 

Without any difficulties one can write similar matrices for rotations about remaining 
axis. 

The biggest advantage of the Stokes vector formalism is the possibility of calcu- 
lating probabilities and cross-sections in a convenient and intuitive way, using scalar 
products and simple matrix operations. The formula 



W=\(l D)(l) 



gives the probability of detecting a particle characterised by the Stokes parameters 
(1 D) in a beam characterised by parameters (IP). In other words, the vector 
(1 D) describes properties of a detector. Of course, a polarisation- insensitive detector 
corresponds to the Stokes vector (1 0); in this case one effectively measures each 



of two orthogonal states and adds the probabilities, leading to the formula: 

W=(l o)( jj). 

Finally, an interaction can be introduced here in a very natural way. When a 
particles experience a polarisation-sensitive interaction, then, in general, their Stokes 
vector is transformed by a 4x4 matrix T, which depends on the interaction type 

W= r U)- (28) 

Thus, the probability of detecting the beam in a state (1 D) after the interaction T 
is given by: 



W 



U^Hp)- 



Full interaction matrices for all processes specific to electrons, positrons and photons 
are available in the Ref. j^j. However, later on in this thesis, reduced matrices (given 
in [TUj ) will be presented and used. The upper left element of a reduced T matrix is 
always equal to one or, precisely speaking, the cross section of an unpolarised beam 
detected by a polarisation-insensitive detector is normalised to unity. 

To summarise the main ideas of this section, we note that: 

• the Stokes parameters describe the polarisation of the beam in a unified way, 

• the interactions can be introduced to the formalism as matrices, 

• (/ P) can be calculated from the positive-energy components of the Dirac equa- 
tion solution, 

• an unpolarised beam can be considered as an ensemble of electrons with spin 
P pointing isotropically in all directions. 

2.3 Passage of electrons through matter 

As it will be described in the sectional the energy, momentum and polarisation of the 
electron from f3~ decay are essential for the i?-correlation determination. However, 
before the electron is detected, it has to cross the whole experimental apparatus and 
in the meantime undergoes numerous interactions with matter. Of course, electron 
energy, momentum and polarisation could have been changed due to these processes 
and hence one has to take them into account in order to retrieve the primary elec- 
tron properties. In the energy region which is in our concern (below 780 keV), the 
dominant processes are ionisation (M0ller scattering) and multiple elastic Coulomb 
scattering. Only a few percent of electrons lose their energy by bremsstrahlung. 

Physics of all these processes is well known and has been already implemented in 
the Geant4, except the depolarisation phenomena. The following section contains 
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some fundamental formulae, which the package is based on (for more details see 
23 ). Additional polarisation effects have been incorporated in the simulation part 
responsible for the electron transportation by means of depolarisation matrices, which 
are given in (TOj and explicitly shown below. 

Last but not least, polarised (3~ decay electrons are Mott-scattered in the analysing 
lead foil, which is actually the main idea of the measurement. Therefore, the last 
section covers the topic of the spin dependent electron scattering on nuclei and its 
analysing power. 



2.3.1 M0ller scattering 

This process occurs when an incident electron is inelastically scattered on an atomic 
electron from the target material. The value of the maximum energy transferable to 
a free electron 1 equals (E incoming — m e c 2 )/2. If the transfered energy is much larger 
than the excitation energy of the material, the atom is ionised and the electron is 
emitted as a so called "<5-ray" . The total cross section per atom for M0ller scattering 



is given by 








,(- FT \ ^^ 


[ 7 2 \2 X ) + x 


1 


27-1, l-x 
7, — m 

7 Z x 


1 — X 


where: 






7 = E/mc 2 , 


(3 2 = 


l-(l/7 2 ), 


x = T cut /(E -mc 2 ), 


y = 


V(7+l), 


r e = the classical radius of the electron. 







T cut = 1 keV is a threshold kinetic energy, below which the process is considered as a 
continuous (in such a case different formulae are used and <5-rays are not simulated). 
The differential cross section, 



da 

de 



27rr 2 Z 



/3 2 ( 7 -l)L 7 
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l) 2 1 



2 7 - 1 

rJ2 



+ 
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27-I 



where e = -wf^ < e < ^, is used to sample the <5-ray energy and direction. With 
respect to the direction of the incoming particle, the <5-ray azimuthal angle <fi is 
generated isotropically and the polar angle 9 is calculated from energy-momentum 
conservation, which corresponds to a situation where target electrons are unpolarised. 
That is how energy and momentum of both, incident and ejected particle is calculated. 
Given these values, one can calculate the depolarisation of the incident electron 
using the formula 12.81 and the reduced depolarisation matrix Tm 



T 



M 



( 1 
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D 
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E 





\o 








F ) 



1 Note that in the M0ller process the scattering and scattered electrons are indistinguishable. 
However, the highest-energy member of the scattered particles is arbitrarily associated with the 
original incoming electron. 
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where: 

C = 2cosfl(27 2 -l)(27 2 -l-7 2 sin 2 fl)/J, 

D = 2-f(2-f 2 -l)sm6cos 2 6/I, 

E = 2cos#(27 2 -l)(27 2 -l-sin 2 fl)/J, 

F = 2[(2 7 2 - l) 2 - (2 7 4 - 1) sin 2 6]/I, 

I = I[(27 2 -l) 2 (4-3sin 2 ^) + (7 2 -l) 2 (sin 4 ^ + 4sin 2 ^)] 

and where 7 is the electron energy (units of mc 2 ) and 9 is the scattering angle, both 
in the centre-of-mass frame (CM). The relations between the CM (not primed) and 
laboratory (primed) quantities 7' and 9' are: 

7' = 2 7 2 - 1, 

. 2 - (Y + 3) sin 2 9' 

cos 9 = - - g — . 

2 + (Y - 1) sin 2 6' 

One can show that the following relation between the corresponding solid angles 
holds: 

[2 + (y-l)sin 2 ^] 2 

The Stokes vector of the incoming particle should be transformed to the coordinate 
system, where the direction of movement is along the z-axis and the xz-plane is the 
plane of the scattering. The polarisation of the final electron is given in the new 
reference frame rotated through the scattering angle 9 about the y axis. 



2.3.2 Multiple scattering 

In addition to inelastic collisions with the atomic electrons, charged particles passing 
through matter also suffer from repeated elastic Coulomb scatterings on nuclei al- 
though with a somewhat smaller probability. Each of them is individually governed 
by the well-known Rutherford formula (when ignoring, for simplicity, spin effects and 
screening) 

ii {d) = ^rX~id^) ^irTf (2 ' 9) 

Moreover, since the nucleus is much more massive than the scattered electron, the 
energy transfer in the process is negligible. 

If the average number of independent scatterings is greater than 20, the prob- 
lem can be treated statistically to obtain a probability distribution for the angle of 
deflection as a function of the particle step length in the traversed material. The 
method governing the multiple scattering (MSC) of charged particles in matter used 
in Geant4 is based on Lewis' theory (see [ijj]) and uses model functions to determine 
the angular and spatial distributions after each step of tracking. The model functions 
have been chosen in such a way as to give the same moments of the distributions as 
the Lewis theory. 



The properties of the MSC process are completely determined by the transport 
mean free paths A&, which are functions of the energy in a given material. The fc-th 
transport mean free path is defined as 



1 f 1 

— = 2wn a / [1 - Pfc(cosx)] 
A& J-i 



-tiT d{cosx) 



where da(x)/dQ is the differential cross section of the single scattering (e.g. the 
Rutherford formula), Pfc(cos x) is the fe-th Legendre polynomial and n a is the number 
of scattering centres per volume. The mean value of the geometrical path length 2 
corresponding to a given true path length 3 is given by 



(z) = A x 



('XI > ( -y 



At the end of the true step length t, the scattering angle is 9, the mean value and 

variation of its cosine are 

t 



{coso) = exp | — — 
Ai 



a 2 = (cos 2 9) - (cos9) 2 = 1 + 2e - e~ 2 \ 

where r = t/Xi and k = t/\2- In addition to this, the square of the mean lateral 
displacement (assuming the the initial momentum is parallel to the z axis) is 



<* 2 + y 2 ) = ^ 



K + 1 K 1 

T + -e -, -re 

K K — 1 K{K — 1) 



KT 



The angular distribution of the scattered particle is sampled according to a model 
function g(u), where u = cosO. The functional form of g is 

g(u) = p[qg x {u) + (1 - q)g a (u)] + (1 - p)g*(u), 

where < p, q < 1, Qi are functions of u, normalised over the range u G [— 1, 1] 



g 1 = C 1 e-« 1 -«\ g2 = C 27K ^, g 3 = Q 



(b-u) c 



^3; 



a,b,c > and Ci are normalisation constants. For more details see the Ref. |2??l 
section 6.2], regarding the purposes of this thesis, it is enough to mention that for 
small scattering angles g\{u) is nearly Gaussian, while for large angles g<i (u) has a 
Rutherford-like tail. 

Finally, the polarisation change in the MSC process is calculated in a very simple 
way (see pHl section 4.3.1]). When the momentum vector of the particle is rotated 
over an angle 0, its component of the polarisation vector in the scattering plane is 
rotated over an angle 6 given by 

s^" 1 



E. 



2 The shortest, straight line distance between the endpoints of a single step. 
3 The path length of an actual particle usually longer than the geometrical path length, since the 
path is random and zigzag. 



where E e is the total energy of the particle in units mc 2 . A simple observation 
indicates that in the higher energy limit, the polarisation of a longitudinally polarised 
beam follows its momentum, while for lower energies, the slower rotation of the 
polarisation vector becomes significant. 

2.3.3 Bremsstrahlung 

For (3~ decay electrons, the radiation of photons in the field of nucleus is much 
less probable process than ionisation or the multiple scattering and affects only few 
percent of electrons, causing minor energy loses and depolarisation. Below 1 MeV the 
theoretical description of bremsstrahlung has to be considered as an approximation, 
the particular model exploited in Geant4 }23t section 7.2] results in up to 15% errors 
for both the cross section and the energy loss. It is based on the tabulated cross 
sections of Seltzer and Berger (Ref. J3|), together with the Bethe-Heitler formula 
(which includes the dielectric suppression of the radiation) and the correction for the 
Landau- Pomeranchuk-Migdal effect 4 . The angular distribution of the emitted photon 
momentum was reported by Tsai J3] and its simplified version is implemented in 
Geant4. 

The depolarisation matrix specific to bremsstrahlung was derived in the Ref. jTTT|. 
section 3.9.2] and is taken as 
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where: 



and: 



I = ( e 2 + e 2 )( 3 + 2 r)-2 ei e 2 (l + 4 M 2 e 2 r), 

E = 4k£Te lU (2£ - 1)//, 

F = 4k£Te 2 u(l - 20//, 

G = 4e 1 e 2 [(l + r)-2u 2 £ 2 r]//, 

H = k 2 [l + %T(£-\) 2 ]/I 



e\ = total energy of the incoming electron (in units mc 2 ), 

e 2 = total energy of the outgoing electron (in units mc 2 ), 

p = electron initial momentum (in units mc), 

k = photon momentum (in units mc), 

u = component of p perpendicular to k, 

u = \u\, 

k = €i — e 2 , photon energy (in units mc 2 ), 



s c 



l 

l+u 2 ' 



4 



The suppression of photon production due to the interference of radiation emitted before and 
after the multiple Coulomb scattering event. 



Both incoming and outgoing polarisation vectors are rotated here into the frame 
defined by the scattering plane (xz) and the direction of the outgoing photon (z-axis). 
Described approach to the depolarisation by bremsstrahlung contains the Coulomb 
and screening effects, introduced in functions: 

T = lnJ-2-/(Z)+^((J/0, 
8 = -*-, 

f(Z) = a?[{l + a 2 )- 1 + 0.20206- 0.0369a 2 + 0.0083a 4 -0.002a 6 ], 

f(Z), where a = aZ, is an approximated form of the Coulomb correction term (see 
jTT]). The function T{8 /£) is tabularised in literature and includes screening effects. 

2.3.4 Mott scattering 

At relativistic energies, the Rutherford cross section ()2.9|) describing the elastic 
Coulomb scattering is modified by spin effects. The resulting Mott cross section 
for the electron, derived from the Dirac equation, may be written as 

do~ M / n \ do~ Rfn ^ ( \ 2 • 2 " 



'^i^in 1 """^ (210) 

however, it gives only the spin-averaged cross section. In order to investigate how the 
particularly polarised electron scatters on a nuclei, one has to introduce the complex 
scattering amplitudes / and g, satisfying the condition 

I(9) = \f\ 2 + \g\ 2 . 

Then, the Sherman function or so-called analysing power S(9) and the spin rotation 
functions T(9) and U(9) are defined as follows: 

S(0) = 1 |^2 + | 3 |f, T(0) = | /p + |J|2, U{6) = |j|2 + | 9 |f- 

As one can see, all observable quantities pertaining to the scattering process can be 
expressed in quadratic terms of / and g. I is the differential cross section for an un- 
polarised beam, S is the asymmetry function which gives the transverse polarisation 
of the scattered electron for an unpolarised beam (or the left-right asymmetry for 
a 100% transversely polarised beam), T and U, finally, describe the rotation of the 
polarisation vector P during the interaction. 

The differential cross section for the polarised electron is dependent not only on 
the angle of deflection 9 but also on the azimuthal angle <fi 

^(e^) = i{9)(i + s(e)P-h) =i(e){i + s(e)p n sin<f>), (2.11) 

where h = n ^g 2 is the unit vector normal to the scattering plane, rii and i%2 are 
unit vectors in the direction of the electron motion, respectively, before and after 
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Figure 2.4: Vectors and angles in Mott scattering (from the Ref.|17j). 

scattering and P n is a component of P along n. After the scattering, the electron 
polarisation is given by 



P 



, _[P n + S (9)] n + T{6)h x[Pxn} + U{Q)[h x P] 

1+P-n 



(2.12) 



The S, T and U functions has been calculated and tabularised for a variety of elements 
and energies by Sherman [T2~] . 
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Chapter 3 
Experiment 



The .R-correlation measurement is being performed at the Paul Scherrer Institute in 
Villigen on the polarised cold neutron beam line of the spallation source SINQ J2j. 
Fig. Kill shows the layout of the line and demonstrates where both the polarisation 
and focusing of the beam take place. 




1. Internal Beam Guide 

2. Beam Shutter 

3. Polarizer 

4. Focusing Guide 

5. Spin Flippers 

E. rvteutron Decay Detection System 



7. Lithium Beam Dump 

8. Electronics 

9. Gas System 

1 D. Radiation Shields 

1 1 . Bo ron ; Lithium Collimator 



Figure 3.1: The spallation source SINQ with the polarised cold neutron beam line. 
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3.1 The main principle 



According to section 12. 1| a measurement of the electron polarisation component, 
which is perpendicular to the plane spanned by the spin of the decaying neutron 
and the electron momentum, will provide an estimation of the f3 decay i?-correlation 
parameter and, if the value was different than zero, would discover time reversal 
symmetry violating processes. The main principle of the measurement is shown in 



Scintillators 




E V 



Wire chambers 



Figure 3.2: The principle of the measurement. The picture shows so called V-Track, 
which is the type of event desired for the transversal polarisation component deter- 
mination. 

Fig. 13.21 An electron from the decay of a neutron with specified polarisation traverse 
the wire chamber detector. Afterwards, it scatters on the analysing foil covered 
with a layer of lead. Note, that the angle of scattering depends on the electron 
polarisation (Eq. 12.11(1 . Deflected electrons traverse back the whole setup, including 
both wire chambers, the second Pb foil and, finally, hit the scintillator, where the 
remaining energy is left. Electron tracks from wire chambers provide possibility 
of reconstruction the emission and Mott scattering angles, which, together with the 
known neutron beam polarisation, is sufficient to obtain the R parameter. In addition, 
the energy measurement from scintillators serves for distinguishing the signal from 
background, generated mostly by high energy electrons coming from neutron captures 
in surrounding. 

The measurements are done for both neutron spin orientations ( "up" and "down" ) , 
which helps to cancel out some systematic effects related to the setup geometry and 
to control better the remaining ones. Moreover, another important feature of this 



experimental configuration is to be used for the same purpose. It is clear that from 
the scattering angles one can acquire information on both transversal components 
of the electron polarisation, because for both of them the experimental technique is 
exactly the same and both can be measured simultaneously. Thus, the iV-correlation 
parameter can be extracted from the data, as well as the R parameter. Since the 
value of N is well known (see Eq. 12.3)1 . the comparison with the measured value will 
serve as an important clue to systematic errors limiting the experimental accuracy. 
Therefore, the behaviour of both electron polarisation components will be studied in 
the simulation. 

3.2 Experimental setup 

The collimated beam of cold neutrons enters the decay chamber filled with helium 
and covered with a special material enriched with 6 Li. Both these elements were 
chosen to reduce the background, because helium has zero cross section for neutron 
capture and is easy to use and 6 Li, which in contrary absorbs neutrons very well, 
does not in effect emit secondary 7 rays, which could convert to electrons in detector 
materials. In the lithium cover, on the sides of the decay chamber there are only 
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Figure 3.3: The helium box with MWPC detectors and hodoscopes. 

two mylar windows separating MWPC (MultiWire Proportional Chamber) detectors 
from decay volume. The MWPC chambers are filled with a gas mixture of helium, 
methylal and isobutan, in proportion optimised for the best efficiency and stable 
operation. The analysing Pb foil is placed behind the next mylar window, in a small 
helium filled box. In this case, helium is meant to protect the foil from oxidising and 
is optimal with respect to energy loss ans multiple scattering. The last elements of 
the setup are two systems of position sensitive scintillators. 
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Chapter 4 
Simulation 



The comprehensive Monte Carlo simulation of the i?-correlation experiment has been 
created by our group since the year 2001 and has already been useful for various 
purposes. Its present version is the sum of individual efforts made by few developers 
and contains all functionality from few recent versions like: 

• EnergyLoss program for energy calibration, 

• Sim used for the investigation of neutron background. 

The first release of the code and its main part, including extended Geant4 classes 
for neutron capture, has been written by J. Pulut. E. Stephan has implemented the 
first version of the Mott scattering subroutine and the code for generation of the 
artificial wire chamber response. 

So far the simulation has been useful for the energy calibration purposes and 
for the optimisation of the experimental setup geometry and materials. Now, with 
the use of new polarised dependent parts, it can be employed for the analysis of 
systematic effects. This chapter describes selected features of the simulation with the 
stress laid on the polarisation transport. 

4.1 The Geant4 package 

Geant4 (GEometry ANd Tracking) is a widely used set of C++ libraries for simu- 
lating the passage of particles through matter, originating in the old Fortran based 
Geant3 version. It includes a complete range of functionality including tracking, 
geometry and physics models and provides tools necessary for implementing a whole 
experiment in a program. With the use of Monte Carlo techniques Geant4 step by 
step calculates a particle energy losses, direction changes, creates secondary particles 
and simulates physical interactions. In general, the total cross section for the vari- 
ous relevant processes is used to select the process that will take place and when it 
is done, appropriate differential cross sections are used to calculate the kinematics. 
However, no electron polarisation sensitivity is included in the code. The details and 
further references can be found in the recent paper of the Geant4 project [22j and 
on the homepage [21__ . 
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Figure 4.1: The geometry implemented in the simulation: a) the collimator, b) the 
proportional chamber, c) the helium box with the analysing lead foil, d) scintillators, 
e) the decay chamber. Directions of axes of the main reference system are shown 
as well. Please note that its origin is located in the middle of the decay chamber, 
between the wire chambers. 



4.2 Geometry 



The whole geometry of the experimental apparatus and the neutron beam imple- 
mented in the simulation is based on its previous version (comprehensively described 
in Ref. [2S]), updated for the new setup. It has been written in a flexible way, there- 
fore, if needed, the user can easily modify materials and sizes of chosen setup elements 
or change the distances between them. All the parameters are available in the macro 
file geom.g4mac presented in appendices and in the source files src/GeometryConstans . ice 
and src/DetectorMaterials. ice . Of course, modification of the source files re- 
quires subsequent recompilation of the program. Fig. 14.11 demonstrates the imple- 
mented geometry and its details. 
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Figure 4.2: The main reference frame {X, Y, Z} and both frames moving with the 
electron. {x p , y p , z p } denotes the particle frame. 

4.3 Polarisation transport 



The electron polarisation a, sampled as described in section 14.3.11 is given relative 
to a right-handed frame moving with the particle. However, the bookkeeping of all 
polarisation changes of a particular electron is performed in a coordinate system 
with the z axis in the direction of motion, in the frame called the particle basis. Both 
frames are linked to each other through a simple rotation, given by a matrix similar to 
that of Eq. 12.71 (see Fig J4.2|) . The full definition of the particle basis {x p , y p , z p } and 
its relation to the laboratory reference system {X, Y, Z} (see Fig. I4.1J1 is as follows: 

• z p is in the direction of particle motion v = (v i, v%, V3), v is normalised 

• x p is parallel to XZ plane 

• Y is parallel to y p z p plane. 

The particle basis definition can be expressed using the formulae 
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4.3.1 Electrons from f3 decay 

The first step of the simulation is the generation of electrons from j3~ decay with 
proper distributions of energy, momentum and polarisation. Although Geant4 pro- 
vides possibility to generate decay of neutrons from the beam, it takes so much 
computational time, that for better efficiency the simulation starts from electrons. 
The neutron beam, however, together with the (3 decay and neutron physics is im- 
plemented in the program and will be helpful in future for identifying main sources 
of background. 

Electrons are uniformly created inside volume of given position and sizes corre- 
sponding to the real shape and spatial placement of the neutron beam. Their kinetic 
energy is sampled using the acceptance-rejection method 1 and an approximated ver- 
sion of Eq. O (see Fig. IP|) 



W(E e )<xp e E e (E„ 



E P 



(4.2) 



Finally, momentum and polarisation are generated from the probability distribution 
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Figure 4.3: Simulated electron energy and momentum spectra (for 1 million events). 



given by 



I¥oc 



AJ)'Ve ^V-Pe ,r-, i^\ ^ &-Pe (J) ' Pe „(J)' (ft X ff) 

1 + A x F + G— P- + Na- (J) + Q — ' + R ^ 

E e E e E e + m E e E e 



(4.3) 

where the polarisation vector a is normalised to unity. It is worth mentioning, that 
values of the decay parameters together with the mean beam polarisation (J) can be 
easily modified by the user (in the file run. mac). In addition to this, user defined 



Also known as the von Neumann method. 
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limits on electron energy and emission angle can be applied which is necessary for 
optimisation purposes. 

Figure l4~4l shows results of a simple test which was performed to check if assumed 
A and G values were reproduced in generated data. Comparing the probability 
distribution 14. HI with the fitted linear function, one can see that 

A s im. = TfT ■ Pi and Gsim. = — 'Pli 

Pe {J} Pe 

where (J) = 0.8973 is the mean neutron beam polarisation and — has to be averaged 
over the (3 decay energy spectrum (what gives (— ) = 1.368). Finally, one obtains: 



A, 



1.368 



•Pi 



-0.118 ±0.003 and G. 



-0.999 ±0.002, 



0.8973 

which is in perfect agreement with the taken values: A = —0.1162 and G = — 1. 
Similar tests have been done for N and R asymmetries in the generated spin of the 
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Figure 4.4: Simulated A and G decay asymmetries with fitted functions 
y = Po{l + pi ■ x) . a is the angle between the electron momentum and polarisation. 
The number of generated electrons equals 1 million. 

electron. Again, comparing the formula 14.31 with Fig. 14.51 after some calculations, 
one can obtain their values straight from the fit. The results are: 

• for assumed R = 0.010 

1.368 



t^sim. 



(J) 



■pi = 0.008 ±0.003, 



for assumed N 



A me 



0.076 



N x . 



— r-pi = 0.072 ±0.002. 
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Figure 4.5: Simulated R and N coefficients. In both cases the fitted function is 
y = p {l+pi-x). 

It is thus confirmed, that j3 decay parameters are reproduced in generated data with 
fair accuracy. 

4.3.2 Depolarisation 
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Figure 4.6: Depolarisation due to the M0ller scattering 

Generated j3 decay electrons traverse the experimental setup, namely the helium 
box and its mylar window, the proportional chamber filled with special gas mixture, 
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Figure 4.7: Depolarisation due to the multiple scattering 
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Figure 4.8: Depolarisation due to bremsstrahlung 



another mylar window and the box containing the Pb foil, also filled with helium, 
where their transversal polarisation is to be analysed by Mott scattering. However, 
before reaching the analysing foil, electrons are partially depolarised through nu- 
merous interactions. This subsection contains some quantitative estimations of the 
polarisation loss based on the simulation results. In order to respect spin effects, 
Geant4 libraries have been extended using formulae from section I2.31 Before the 
actual results have been acquired, an especial effort had been made to check the 
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reliability of the program and to verify the outcome with Ref. [TU1 section 4.3]. 

The figures below (j4.6l - l4~Hj) show the depolarisation effects due to three imple- 
mented processes, separately (Geant4 gives possibility to "turn off" chosen physical 
processes). One million fully transversely polarised electrons have been shot from the 
point exactly in the middle of the helium box with momentum vector perpendicular 
to the Pb foil surface and the kinetic energy set to 700 keV. The polarisation vector 
in the particle frame was P = (1,0,0) in order to investigate behaviour of the polar- 
isation component strictly related to the i?-correlation. The "final polarisation" on 
the pictures is the P x component at the surface of the analysing foil, after the trans- 
portation through the whole setup. As it was expected, ionisation is the dominant 
process and the least change is caused by bremsstrahlung. 

4.3.3 Mott scattering in the analyzing Pb foil 

After crossing the proportional chamber, a partially depolarised electron can hit the 
analyzing foil covered with a layer of Pb, where it is scattered (as described in section 
I2.3.4J) . The formula 12.111 which governs this process, can be rewritten in a form 

^(0A) = HO) [1 + S(6) (P y cos <j> - P x sm<f))} , (4.4) 

which makes the dependence on the transversal polarisation component more obvious. 
As it was already mentioned, Geant4 does not contain any electron polarisation 
dependent processes, therefore to include the Mott scattering some extensions of the 
code have been done. The most important steps of the implemented routine are: 

1. Functions / and 5* are loaded from the file (see Fig. I4.9|) . 

2. As soon as an electron reaches the surface of the analyzing foil, the Geant4 
tracking engine is stopped. 

3. An exact point of the scattering inside the foil is generated from the uniform 
distribution. 

4. The energy loss between the foil border and the scattering point is calculated 
using the database [21] and subtracted from the electron energy. 

5. The new particle momentum is generated from the formula 14.41 If needed, the 
values of functions / and 5* are linearly interpolated. 

6. The energy loss between the scattering point and the foil border is calculated 
and subtracted. 

7. The new particle momentum and energy are passed to the Geant4 tracking 
engine, which is started again. 

And here appears an important issue. The formula 14.41 is valid for the scattering 
on a single nuclei, what corresponds to the exact point of the scattering which is 
sampled in the third step of the routine. However, the electron in the foil can be 
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also scattered two or more times, what cannot be neglected, even though for the 
foil thickness used in the experiment the single scattering is still more probable (the 
thicker is the Pb layer, the more probable is double or multiple scattering). The 
solution of this problem are, calculated for a given foil thickness, effective functions 
I e ff9) and S e a{6) which already contain multiple scattering effects and ought to be 
used instead of the normal cross section and analysing power. Both effective functions 
have been calculated by E. Stephan J4J with a dedicated Monte Carlo simulation based 
on Ref. JH] and are presented in Fig. 
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Figure 4.9: The differential cross section and the effective Sherman function used 
in the program. Both are results of an independent simulation performed by 
E. Stephan gj. 

In order to check reliability of the Mott scattering routine a simple test have been 
prepared. Scattering of one million fully transversely polarized electrons with the 
polarisation vector 



-* — \*X7 *-yi *z 



vVV5' 



= (0.447,0.894,0) 



and kinetic energy 700 keV have been simulated, resulting in angular distributions 
shown on Fig. 14.101 The scattering angle 9 S here, and for all other simulations, was 
restricted to the region 120 - 160 degrees, because for lower angles reflected electrons 
do not have a chance to hit the detector and for larger angles the cross section goes 
down rapidly. The cross section 14.41 averaged over the 9 S angle can be written as 

^(0 S ) H^> [1 + (S> sin (0 S + 5)] , (4.5) 

therefore full information about the transversal polarisation component is hidden in 
the phase shift 5 and can be obtained from the fit (see Fig. 14.101) . The results 



P x = - cos 5 = 0.454 ± 0.006 P v = sin 5 = 0.891 ± 0.003, 
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are in a very good agreement with assumed values, what provides a satisfactory 
crosscheck. 
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Figure 4.10: Test of the Mott scattering routine: distributions of the scattering angles. 
The function fitted to the S distribution is y = (I) [1 + (S) sin(0 s + 5)]. 



One should be aware of several systematic errors that might be introduced here 
and could lower the accuracy of the simulation, especially when it will be used for 
prediction of the real measurement results: 

• the effective functions I & a and S & a are model dependent and can contain some 
systematics caused by the generation method ([T%]). 

• because of the unavailability of proper data for 82Pb, the analyzing power S for 
soHg has been used as the input for the S e a calculation, 

• small additional deviations are caused by the linear interpolation of I e « and 
S e a values between the points of given energy and scattering angle, 

• so far the electron depolarisation inside the Pb-foil, directly before the scattering 
has not been implemented. 



Chapter 5 
Systematic effects 



After numerous careful tests of the program, described in the previous chapter, it 
finally could have been employed in realistic simulation of the whole experiment 
and its conditions. The results presented below required two days of calculations 
on a PC machine with a fast CPU. It was enough to reach the expected statistics 
of the planned experiment, namely, over one million events (V- Tracks). Due to the 
significant differences in cross sections and analysing powers, the data analysis was 
performed separately for three electron energy intervals: 

• 200 keV - 400 keV, 

• 400 keV - 600 keV, 

• 600 keV - 800 keV. 

Electrons with energies below 200 keV were not generated, since the threshold of 
our apparatus is expected to be around this value. For optimisation reasons, after 
dedicated tests, some cuts have been also applied on electron generation angles: 

• the polar angle: 9p E (15°, 165°), 

• the azimuthal angle: (pp E (20°, 160°). 

Moreover, the Mott scattering angle 9 S was restricted to a subset (120°, 160°). The 
motivation is clear, on the one hand side, there is no chance to detect in the opposite 
chamber an electron scattered with 9 S < 120°, on the other hand for the scattering 
angles over 160° electron tracks before and after scattering are not distinguished by 
the wire chamber. 



5.1 Data analysis 

First, it is necessary to describe the general concept of the data analysis. For each 
Mott scattering event (V- Track v ) the simulation saves to a file the data specified 
below: 



*An electron scattered from the foil which hits the scintillator on the opposite side of the beam 
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• initial momentum direction, 

• momentum at the Pb foil surface, before the scattering, 

• momentum after scattering, 

• initial polarisation, 

• polarisation at the Pb foil surface, before the scattering, 

• the initial electron position, 

• the scattering point position, 

• final position at the scintillator border, 

• initial kinetic energy, 

• final kinetic energy, 

• random deviation from the final kinetic energy (to simulate the effect of the 
finite energy resolution), 

• artificially generated response of electrodes from the wire chambers. 

The raw data files, containing values in the ASCII format, are later converted with the 
extended version of the NPRUN0 analysis program to binary *.root files. During 
this process a few important parts of the actual analysis are done, including the 
reconstruction of V- Tracks from the artificial wire chamber data. Afterwards one 
can work directly on the binary files, using the ROOT |2Zj environment and macros, 
created especially for that purpose. At this point, the analysis is done separately 
for momenta and positions reconstructed from artificial tracks and separately for the 
"real" values, taken straight from the simulation. 

To obtain values of R and N coefficients from the simulated data, one has to know 
the intensity distribution of the azimuthal scattering angle (f> s , which can be used to 
calculate mean values of the transversal polarisation components, just as it was shown 
in section 14.3.31 Both scattering angles can be calculated from momentum directions 
before and after scattering, transformed to the incoming particle frame (see section 
I4.3J1 . The Fig. 15.11 shows the distribution of scattering angles for over 200000 V- 
Tracks, generated from the beam polarised in the "up" direction. Its shape is mostly 
a consequence of acceptance of the the experimental setup and beam geometry and in 
this form cannot be used for determination of the electron spin direction. To achieve 
a distribution independent on the geometry, just like in the real experiment, one has 
to flip the beam polarisation and produce the same amount of data. In this case, 
the influence of the geometry is exactly the same, while the effects caused by the 
polarisation contribute to the <p s distribution with the opposite sign (from equations 
14.31 and 14. 4jl . After subtracting both distributions and normalising the result to the 
total number of events, one gets the final, geometry independent intensity distribution 
of Mott scattering angles for V- Tracks. The mean values of polarisation components 
can be now extracted from the fit. 
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Figure 5.1: The comparison of the distribution of real scattering angles <p s (pictures 
on the left) with the reconstructed ones (on the right). 

5.2 Reconstruction of V- Tracks 



Probably the most challenging problem in the analysis of data from the real ex- 
periment is the reconstruction of electron tracks based on signals from the MWPC 
chambers. To avoid depolarisation and multiple scattering effects, the chambers are 
relatively thin (around 10 cm) and contain only 5 planes of electrodes. Each plane 
consists of one layer of anodes (horizontal wires) and two layers of cathodes (vertical 
wires) . 

The reconstruction algorithm, implemented in the NPRun program is in its ad- 
vanced development phase and still requires some tests, while the simulation pro- 
vides the only possibility to compare "real" tracks with the reconstructed ones and 
to check the reconstruction efficiency Therefore, some qualitative comparisons have 
been done. As we can see in Fig. 15.11 the reconstruction works fine, with very high 
efficiency, however, for the angles <p s = 0, ±~, ±7r there is a large difference. The miss- 
ing reconstruction efficiency in that specific geometry is due to MWPC wire layout 
feature. The reconstruction of V- Tracks has to be done separately for two projec- 
tions (for anodes and cathodes). The problem arises for V- Tracks which lay either in 
the cathodes or in the anodes plane. Since they can be fully reconstructed in only 
one projection, ambiguities can appear. The scattering angle values <p s = 0, ±|, ±7r 
correspond exactly to these events. The version of the algorithm tested in this case 
reconstructs as much of these type of V- Tracks as possible which results in many 
misidentified tracks. At the final stage of the analysis the effect cancels out (on the 
cost of lower statistics) and the ultimate result is satisfactory (Fig. I5.2J) . It should be 
added here, that in reality the efficiency of the reconstruction algorithm is much lower, 
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Figure 5.2: The comparison of the real scattering angles (picture on the left) with 
the reconstructed ones (on the right) after the subtraction of data corresponding to 
opposite beam polarisations. 

due to the noisy signal from the chambers caused by various sources of background. 
The algorithm for the artificial chamber signal generation includes the optional noise 
generation, but in this case it has not been used. 



5.3 False R-correlation 



Although, the false contribution to the R coefficient caused by the neutron decay 
asymmetry was one of the expected systematic effects (a similar effect appeared in 
the R measurement for 8 Li [IE]), its nature have been understood and the magnitude 
have been estimated very recently, by means of this simulation. Figures 15.31 and 
15.41 explain the source of the effect. Let's consider the situation when the beam 
is polarised in the "up" direction. For simplicity, the dependency on the electron 
polarisation is not taken into account. Because of the neutron decay asymmetry 
(nonzero A coefficient), more electrons are emitted "down", in the direction opposite 
to the beam polarisation. In the case of the electrons emitted "up" , the scattering 
angle value (fi s = — | is less probable than |, since only the backs cattered particles 
are accepted by the simulation. Knowing that in the program 9 S G (120°, 160°), 
one can see that the effect starts to appear for the angle between the incoming 
electron momentum and the scattering foil plane a < 60°. In contrary to this, for 
the much larger number of particles emitted in the direction "down", the effect is 
opposite and <p s = — ~ is more probable than |. Of course, if the beam polarisation 
is flipped, the effect is reversed, therefore it cannot be canceled by subtraction of S 
distributions measured for both beam polarisations. It results in the final sine-like 
angular distribution, with one maximum at <p s = — | and one minimum at (f) s = ~ 
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Figure 5.3: The neutron beam as seen from the position at the analysing Pb foil. 
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Figure 5.4: The principle of the systematic effect. For electrons emitted in the 
direction "up" the scattering angle S = — | is less probable than the angle |. 
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Figure 5.5: The false R effect for the initial electron energy Ek £ (200, 782)keV. N u 
and N d denote the count rate for the "up" and "down" beam polarisation, respec- 
tively. The function y = Pq- sin(0 s + 5) has been fitted to data points. 



(see Fig. I5.5J) . As it was specified before, it is exactly the shape that one would 
expect for the nonzero R coefficient. Therefore, it is necessary to examine this effect 
and estimate the induced false contribution to the measured R coefficient. 

A dedicated simulation have been performed for that purpose. In order to separate 
the effects induced by the f3 decay asymmetry, all the coefficients except A have been 
set to zero. The program has produced almost 2 million V-Tracks, half of them with 
the beam polarisation "up" and half with the opposite. The total result averaged over 
the electron kinetic energy Ek > 200 keV and over all generation angles is presented 
in the Fig. I5.51 However, even more interesting are similar pictures for three distinct 
electron energy intervals (Fig. I5.6|) . The amplitudes p$ taken from the fit can be 
now used to calculate the false contribution to the average value of the transversal 
electron polarisation component. The intensity distributions are given by: 



N d 



(I) [1 + (S) sin(0 s + 5)] = (I) [1 + (S) (P y COS0, - P x sin <j>,)], 



(I) [1 + (S) 



-Py COS 



+ P x sin(j) s )}. 



thus 



N u - N d 2 (I) (S) {P y cos <p s - P x sin . 



N u + N d 



2(1) 



(S)smU s + 5). 
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Figure 5.6: The false R effect for three electron energy intervals. The function y 

Po ■ sin(0 s + 5) has been fitted to data points. 
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The average Sherman function values for specified energy ranges are: 



E k [keV] 



(S) 



200 - 400 



-0.15 



400 - 600 



-0.20 



600 - 800 



-0.23 



and both transversal polarisation components: 

cos 5 



The results are: 



P, 



E k [keV] 



(S) 



•Po, 



P, 



sin 5 



■Po- 



P.,: 
P„ 



200 - 400 



0.14 ±0.01 
0.00 ±0.01 



400 - 600 



0.13 ±0.01 
0.00 ±0.01 



600 - 800 



0.11 ±0.02 
0.00 ±0.02 



Note, that the systematic error of the average S has not been taken into account in 
the error estimation. An obvious conclusion from the above simple analysis is that 
the magnitude of the investigated systematic effect does not depend on the electron 
energy. 

It should, however, depend on the electron emission angle as it was explained 
before. And that was the motivation of the last test that has been done. The false 
R effect has been extracted from the data separately for different values of the polar 
angle between the normal to the scattering foil (the x axis in the main reference 
system) and the initial electron momentum. In this case, events with all possible 
initial electron energies have been treated together and the average value of analysing 
power for the whole energy range 200 - 782 keV was assumed to be (S) = —0.18. 
The result is shown in the Fig. 15 .7[ one can see a clear dependency on the emission 
angle, what confirms that the nature of the systematic effect is well understood. For 
the emission angle below 50° its magnitude increases linearly (a similar situation like 
in Ref.|15j). later it starts to decrease. It can be explained by the fact that for the 
larger emission angles electrons may scatter on the MWPC frame, therefore their 
final angle at the foil can be dramatically changed. Accurate estimation of the false 
contribution to the R coefficient is essential for the final systematic uncertainty of 
the whole measurement. 



O.d. f/\lj^tlj n-UUfi,ltfiLAliUi\ 



4/ 



C 

o 

c 
o 

Q. 

E 
o 
o 

c 

'q. 

"5 
w 

i- 

(U 

> 
c 

CO 



0.2 










































0.15 


























0.1 
























0.05 
-0 
0.05 
-0.1 
0.15 








































































-, , , , 


i , , , 


, , , , 


, , , , 


, , , , 


, , , , 


, , , , 


, , , 


( 


) 


1 


2 


3 


4 


5 


6 





70 



The electron emission angle 



Figure 5.7: The false R effect as a function of the electron emission angle. 



45 Lyn/\ril^n 0. DIDliiiVi/lIiLy il/l^il/L/ID 



Chapter 6 

Summary and conclusions 



The main purpose of this thesis was the development of the comprehensive Monte 
Carlo simulation and the analysis of systematic uncertainties in the planned R- 
correlation coefficient measurement. The biggest effort was required to create and 
test the simulation program, based on the Geant4 package and to develop appro- 
priate data analysis tools. In the first step, all previously written parts of software 
have been checked, debugged and unified in one application providing the full func- 
tionality. Furthermore, the new construction of the experimental apparatus had to 
be included replacing the old implementation. The most challenging part of the de- 
velopment phase was, however, the necessity of creating the code for the simulation 
of polarised electron processes and for the polarisation transport. So far, such an 
option has not been provided by the authors of Geant4 package. 

Finally, the program has been used to search for systematic effects in the mea- 
surement. The first effect associated with the neutron /3 decay asymmetry has been 
already recognized, explained and estimated. Due to the false contribution to the 
R coefficient induced by this phenomena, it is essential to control it as well as it is 
possible. The Monte Carlo simulation provides the robust and powerful tool, which 
will serve for that purpose. Another important issue is the problem of the electron 
track reconstruction, which might introduce new systematic effects of different type. 
With artificial data, these contributions can be better controlled. 

The search for other effects is in progress and definitely the program is going to 
be extensively used in future for the further measurement improvement and the data 
analysis. 
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Appendix A 



The macro run. mac 



/gun/PositionRandom on 
/gun/PositionRandomType linear 

/gun/EnergyRandom on 
/gun/EnergyRandomType raw 
/gun/Decay on 
/gun/particle e- 

/gun/MomentumRandom on 
/gun/AngleMinus 70 degree 
/gun/AnglePlus 70 degree 

# Neutron beam properties 
/gun/XBeam 45. mm 
/gun/YBeam 145 . mm 
/gun/XDiv 0.8 degree 
/gun/YuDiv . 9 degree 
/gun/YdDiv . 65 degree 
/gun/range 1 . 1 

/gun/neutron_polarization 0.8973 
/gun/NeutronSpinRandom on 

# Neutron beta decay parameters 
/gun/A 0.0 

/gun/R 0.0 
/gun/N 0.0 

/chamber/threshold_v .4 
/chamber/threshold_h .8 
/chamber/sjit_v 10. 
/chamber/sjit_h 10. 
/chamber/rnoise_level_v 0. 
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/chamber/rnoise_level_h 0. 
/chamber/ef fiplane 0.95 



/run/beamOn 1 
#/run/beamOn 50000000 



Appendix B 

The macro geom.g4mac 



#Titan Foil 

#/mydet/boron_box/ti_f oil_thickness 5.393 micrometer 

# full thickness of the Ti foil (for energy loss simulation) 
/mydet/boron_box/ti_f oil_thickness 0.0 micrometer 

# Boron box 

/mydet/boron_box/bigger_x 700. mm # boron box width (with at the end of the box) 
/mydet/boron_box/smaller_x 214. mm # boron box width (the central part) 
/mydet/boron_box/height 625. mm # height of the b.b. 
/mydet/boron_box/length 2533.5 mm # total b.b. length 

/mydet/boron_box/begin_length 360. mm # length of the beginning of the b.b. 
/mydet/boron_box/end_length 1278.5 mm # length of the ending of the b.b. 
/mydet/boron_box/border 1.25 mm # total b.b. border thickness 
/mydet/boron_box/mid_border 1.2 mm # thickness of the inner layer of the border 

/mydet/boron_box/window/height 500. mm # height of b.b. windows 
/mydet/boron_box/window/width 500. mm # width of b.b. windows 

# window placement : 
/mydet/boron_box/window/y_pos 0. cm # - height 
# (relative to the geometrical center of the b.b.) 
/mydet/boron_box/window/z_pos 810. mm # - z-axis position (along the beam) 

# Dump 

/mydet/boron_box/dump/width 650. mm # beam dump sizes 
/mydet/boron_box/dump/height 622.5 mm # 

/mydet/boron_box/dump/thickness 3.6 mm # and its placement 
#/mydet/boron_box/dump/position 2530. mm # along the beam line 

# MWPC 

/mydet/mwpc/height 625. mm # MWPC sizes: 
/mydet/mwpc/length 625. mm # 
/mydet/mwpc/thickness 91. mm # 
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/mydet/mwpc/zposition_left 810. mm # rel . to the beam: left chamber 

/mydet/mwpc/yposition_left 0. cm # 

/mydet/mwpc/zposition_right 810. mm # right chamber 

/mydet/mwpc/yposition_right 0. cm # 

# Wire chamber 

/mydet/mwpc/inside_height 500. mm # chamber sizes (inside the MWPC) 
/mydet/mwpc/inside_length 500. mm # 

/mydet/mwpc/anode_spacing 5. mm # distance between anodes 

/mydet/mwpc/dist_anode_cathode 4. mm # distance between anode and cathode planes 
/mydet/mwpc/cathode_spacing 2.5 mm # distance between cathodes 
/mydet/mwpc/anode_plane_spacing 16. mm # distance between anode planes 
/mydet/mwpc/anode_diameter 25. micrometer # wire diameter (anodes) 
/mydet/mwpc/cathode_diameter 25. micrometer # wire diameter (cathodes) 
/mydet/mwpc/dist_entr_anode 12. mm # distance between the chamber 

# entrance and the nearest anode 
/mydet/mwpc/anode_planes # total planes number 

/mydet/mwpc/material/isobutane 10. # gas mixture: -percent of isobutane 
/mydet/mwpc/material/methylal 6. # -percent of methylal 

# Scintillators 

/mydet/scintillator/height 60. cm # scintillator sizes 
/mydet/scintillator/width 10. cm # 
/mydet/scintillator/thickness 10. mm # 

/mydet/scintillator/spacing 1 . mm # scintillator spacing 
/mydet/scintillator/number 6 # number of scintilators 

# Golden (lead ?) foil 

/mydet/f oil/width 60. cm # foil sizes 

/mydet/f oil/height 60. cm # 

/mydet/f oil/foil in # "in" or "out" (with and without the foil) 

# Distances 

#/mydet/to_left_mwpc 0. mm # between b.b. surface and the left MWPC window 
/mydet/to_right_mwpc 0. mm # between b.b. surface and the right MWPC window 
/mydet/to_scintillator 6.6 cm # between MWPC's and hodoscopes 
/mydet/to_golden_f oil 3.2 cm # betweem MWPC's and the lead foil 

# Nose (Collimator) only rough values 

/mydet/nose/bars_number 6 # number of barriers inside the collimator 
/mydet/nose/to_barrier 50. mm # distance between the nose ending and 

# the first barrier 

/mydet/nose/border 11.2 mm # total thickness of the "nose" border 
/mydet/nose/mid_border 1.2 mm # thicknes of the inner layer of the border 
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